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The atom-dimer and dimer-dimer scattering lengths are analytically calculated in an expansion around four 
spatial dimensions for fermions with a large 2-body scattering length a. We find the atom-dimer scattering 
length dad /a = j — |e + 0(e 2 ) and the dimer-dimer scattering length a^ja = £ — 0.344e 2 + 0(e 3 ), where 
e = 4 — d and d is the number of spatial dimensions. These ratios a ac / /a = 1.11 and ajj ja = 0.656 at £ = 1 
are to be compared with the non-perturbatively calculated numerical results 1.18 and 0.6 respectively. The 
neutron-deuteron scattering length in the quartet channel a„o ~ 4.78 fm is in reasonable agreement with the 
experimental value = 6.35 ±0.02 fm, considering 2-body effective range corrections ro/a ~ 0.4 were 
excluded. The deuteron-deuteron scattering length cidd ~ 3.15 fm in the spin-2 channel. 



Systems with a large 2-body scattering length a are com- 
mon in nuclear and atomic physics. Recent experiments 
near the Feshbach resonance have added new systems to this 
list Ql|2>|3l0,|5|. Near the unitarity limit a — > ±°o, the leading 
behavior is expected to be universal on dimensional ground. 
However, theoretical calculations in this interesting regime of 
large 2-body scattering length are difficult - typically involv- 
ing non-perturbative resummation of diagrams or non-trivial 
numerical calculations. In Ref. |6], Nussinov and Nussinov 
pointed out that in d = 4 spatial dimensions one obtains a non- 
interacting system of dimers. Recently Nishida and Son have 
developed an e = 4 — d expansion around four spatial dimen- 
sions l7j], allowing them to perturbatively calculate thermody- 
namic properties for a dilute Fermi gas at the unitarity limit 
a — > ±°o. 

In this paper, we consider the calculation of the atom-dimer 
and dimer-dimer scattering amplitude at large but finite 2- 
body 5-wave scattering length a such that the interaction is 
completely determined by 2-body physics. Besides providing 
a perturbative method for few-body calculations in the vac- 
uum, this work would allow study of many-body properties 
at finite a such as the Bose-Einstein condensation (BEC) in 
atomic systems near the Feshbach resonance. The relevant 
theory for large scattering length a is described by the follow- 
ing Lagrangian density: 
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where \\i are spin-^ fermionic fields and the pauli matrix c?2 
projects the interaction onto the singlet 5-wave. An equivalent 
form for the Lagrangian density after a Hubbard-Stratonovich 
transformation is 
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In d = 3 spatial dimensions, it is possible to establish a 
systematic perturbation in R/a where the range of the in- 
teraction R <C a |H Q3 fllll . This has been quite suc- 
cessful in the 2 and 3-nucleon systems, see for example 
Refs. 0, where the 2-body scattering length is 
large. However, it involves resumming certain diagrams to 
all order in perturbation. In Fig. ^ tne atom-atom scatter- 
ing amplitude is 0(a) and requires resumming an infinite set 
of diagrams all of which are leading order in the perturba- 
tion 13 Q3 Hill . The leading order atom-dimer amplitude 
is again 0(a) and given by an infinite set of diagrams, see 
Ref. HI. The dimer-dimer amplitude is also 0(a), and sim- 
ilarly requires resumming the indicated diagrams at leading 
order |15]. The atom-dimer amplitude for spin-i fermions 
has been non-perturbatively calculated model-independently 
in effective field theory llflflZfloll . For dimer-dimer scatter- 
ing only a partial resummation of the diagrams has been done 
in Ref. (Tj. 
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FIG. 1: The leading order scattering amplitudes for atoms and 
dimers in d = 3 spatial dimensions. Single lines: atoms/fermions, 
double lines: fully dressed dimers/bosons. 
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Working in d = 4 spatial dimensions will avoid the non- 
perturbative treatment of the dimer mentioned above j3l- We 
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will formulate an e expansion for atom-dimer and dimer- 
dimer scattering in d = 4 — e dimensions |7[ ■ First we define 
the coupling co in arbitrary d spatial dimensions using the pole 
in the 2-body scattering amplitude, Fig.^ 
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located at the 2-body binding energy B=l/(a 2 M): 
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and then write the Lagrangian density in terms of the kinetic 
and interaction piece as 
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The original Lagrangian density in Eq. Q had only one 
coupling cq that was determined by the 2-body scattering 
length a, the only parameter in the theory. In Eq. 0, after the 
Hubbard-Stratonovich transformation, we have an additional 
coupling g. However, its value is arbitrary and only the ratio 
of the yukawa coupling squared g 2 and the "mass" term g 2 /co, 
i.e. co, contributes after the dimer wavefunction renormaliza- 
tion. In d = 3 dimensions, the wavefunction renormalization 
for the dimer Z</ ~ I /a II ill . In d = 4 — e dimensions, where 
we expect a non-interacting theory of dimers 0| , the coupling 
g is chosen such that using the Lagrangian density in Eq. l|6) 
gives Zd ~ 1 + 0(e)- This then requires the dimer self-energy 
fL(pQ,p) in Fig.|2]to vanish on-shell and only contribute to Zd 
through a derivative term at 0(e): 
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where "• • • " indicate terms higher order in the e expansion. 

In e = 4 — d dimensions, for Eq. to hold, the coupling g 
has to be 
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FIG. 2: Dimer wavefunction renormalization. Solid lines: 
atoms/fermions, dashed lines: dimers/bosons. Unlike in d = 3 di- 
mensions, the dimer propagators are not dressed since the dimer- 
fermion interaction is perturbative around d = 4. The "X" repre- 
sents insertion of ig 2 /cq and the "square" kinetic energy —i[po + 
B — p 2 / (2M)] from £/ in Eq. jfj, on the dimer propagator. 



where only the factor £ e is uniquely determined from the 
dimer wavefunction renormalization. With this choice of the 
coupling g, it is straightforward to estimate the relative sizes 
of diagrams. For elastic scattering, diagrams involve even 
powers of g and we associate a power of e for every factor 
of g 2 . In effect, the e expansion is a loop expansion. 
For the wavefunction renormalization Z e i, we get: 
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where Je ~ 0.57722 is the Euler-Mascheroni constant. 

The atom-atom scattering amplitude is given by the dia- 
grams in Fig. |3] We find at threshold: 
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where we divided the amplitude i7 aa by factors a 3 ~ e to look 
at dimensionless ratios for convenience. Unlike in d = 3 di- 
mensions, the atom-atom scattering amplitude i1 aa scales as 
0(a 2 ) in d = 4 dimensions. 
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FIG. 3: Atom-atom scattering amplitude 7 aa in the e expansion. Di- 
agrams up to the next-to-leading order are shown. Same notation as 
in Fig. [2] 

A non-trivial check of the e expansion with finite a would 
be the atom-dimer scattering length. Like the atom-atom scat- 
tering amplitude, the atom-dimer amplitude i1 a( i also starts at 
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0(e), see Fig.|4] At 0(e 2 ), there is a contribution from a 1- 
loop diagram shown in Fig.0] in addition to the wavefunction 
renormalization factor Z c / calculated in Eq. l|9}. The contri- 
bution from the tree diagram (including a factor of a E ~ 3 ) in 
Fig. His: 
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Together with the contribution from the 1-loop diagram, we 
get at next-to-leading order 
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To compare with numerical results calculated in 3 spatial di- 
mensions, we consider the dimensionless ratio 
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In d = 3 spatial dimensions, 7 a d /7 a a = 3a a d / (4a). Numerical 
evaluation by solving the 3 -body Faddeev equation flol fl2l 
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gives 7 a d/7 aa ~ 0.885 (a a d m 1.18a). The e expansion result 
is perturbatively close to the non-perturbatively calculated nu- 
merical value, and it is a tremendous simplification over the 
non-perturbative calculation in d = 3 dimensions. 
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FIG. 4: Atom-dimer scattering amplitude 7 ac i in the e expansion. 
Diagrams up to the next-to-leading order are shown. Same notation 
as in Fig.0 

The discussion can be generalized to certain spin chan- 
nels in nuclear physics where one has to introduce isospin 
degrees of freedom to accommodate both the proton and the 
neutron with similar properties under the strong nuclear inter- 
action. Specifically, the dimer in nuclear physics context is the 
deuteron, a bound state of a proton and a neutron, in the spin- 1 
channel. The spin degrees of freedom in the dimer calculation 
can be thought of as the isospin of nuclear physics, and then 
the dimer calculation can be applied to deuteron scattering in 
channels where all the proton and neutron spins are pointing 
in the same direction. Thus the atom-dimer amplitude is the 



same as the neutron-deuteron scattering in the 5-wave quartet 
channel with angular momentum J = 3/2. 

The deuteron binding energy of B — 2.2246 MeV corre- 
sponds to a = 1 / \JMB = 4.32 fm, ignoring the effective range 
r{) = 1.75 fm 1131 . The standard nuclear physics definition of 
the scattering length a ~ a + ro/2 sa 5.2 fm. From the atom- 
dimer calculation, the neutron-deuteron quartet channel scat- 
tering length is a n o ~ 1.11a w 4.78 fm. This is in reason- 
able agreement with the experimental value = 6.35 + 0.02 
fm lllil considering we neglected ro/a ~ 0.4 effects. For 
the deuteron, the effective range ro corrections are significant 
even though formally higher order in the perturbation. The so- 
called Z-parameterization was developed in d = 3 dimensions 
to precisely account for these effects 12011 . 

Finally, the dimer-dimer scattering length is calculated in 
the e expansion. The leading contribution to the dimer-dimer 
scattering is 0(e 2 ) from the 1-loop diagram shown in Fig. [5] 
We calculate the dimer-dimer scattering amplitude to next-to- 
leading order where it receives a 2-loop contribution. The 1- 
loop contribution is: 
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The 2-loop contribution gives 
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Thus, for the next-to-leading order dimer-dimer amplitude, we 
get the analytic result: 
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7dd/7aa = add /2a in a" = 3 spatial dimensions and numeri- 
cal calculation gives add = 0.6a ll2lll . a value perturbatively 
reproduced in the epsilon expansion. From the dimer-dimer 
calculation, the deuteron-deuteron scattering length in chan- 
nels with all the proton and neutron spins in the same direction 
is aoD ~ 3.15 fm in the e expansion, ignoring 2-body effective 
range ro effects. 

Systems with unnaturally large 2-body scattering lengths 
are common in nature and typically necessitate non- 
perturbative calculations that are computationally difficult and 
expensive. Study and development of perturbative techniques 
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FIG. 5: Dimer-dimer scattering amplitude 7^ in the £ expansion. 
Diagrams up to the next-to-leading order are shown. Same notation 
as in Fig.|2| 

are always useful. Many times perturbative calculations can 
be done analytically that allow insights otherwise not possi- 
ble in a numerical calculation. In this paper, we considered a 
fermi system with a large but finite 2-body scattering length a 
in an expansion around d = 4 spatial dimensions. The effect of 
the 2-body scattering length was systematically incorporated 
and the atom-dimer and dimer-dimer scattering length was 
calculated to next-to-leading order in perturbation. The results 
seem to converge to the non-perturbatively calculated numer- 
ical solutions at this order of the calculation, and the e expan- 
sion should be an useful tool for studying strongly interacting 
systems with large 2-body scattering length a. For example, 



the current work would be important for research in the BEC- 
BCS crossover region in atomic systems as one varies the scat- 
tering length a across the Feshbach resonance. As shown, 
the dimer calculation can be generalized to deuterons in nu- 
clear physics in certain spin channels. However, it is not clear 
from this calculation how the Efimov effect l22ll in bosons and 
fermions with isospin (more than one fermion species as in 
the proton-neutron system) associated with non-perturbative 
renormalization |23] emerges in the e expansion. It could be 
that the Efimov effect cannot be described by an analytic con- 
tinuation from d = 4 spatial dimensions. It is known that the 
Efimov effect occur only in spatial dimensions 2.3 < d < 3.8, 
with d = 3 being the only integer dimension l24ll . 
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